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i.  INTRODUCTION 


Target  scintillation  is  the  variation  of  angular  location  of  the  apparent 
source  of  the  echo  signal  and  is  dependent  upon  relative  phase  and  amplitude 
of  the  component  echo  signals  and  their  apparent  angular  locations.  The 
motion  of  the  target  causes  the  apparent  source  to  wander  back  and  forth 
about  the  physical  center  of  the  target.  This  wander  is  called  angle  noise  or 
scintillation  [Ref.  1 ] . Since  the  noise  is  random,  it  is  usually  calculated 
using  experimental  observation  data.  Few  articles  have  appeared  to  propose 
models  of  the  phenomenon. 

Recently,  we  implemented  a simulation  of  target  scintillation.  The 
objective  of  this  report  is  to  document  this  implementation  and  to  define  its 
theoretical  basis.  These  objectives  are  realized  by  answering  the  following 
questions: 

1.  What  algorithm  was  actually  implemented. 

2.  What  are  the  statistical  characteristics  of  simulated  scintillation? 

3.  What  are  the  theoretical  statistical  characteristics  of  our  model? 

2.  ALGORITHM 

The  scintillation  implementation  is  based  on  the  model  illustrated  in 
Figure  1 . Random  noise  is  generated  by  the  algorithm  derived  in  Appendix  A 
The  noise  is  passed  through  a filter  to  simulate  the  correlation  of  scintillation. 
T’he  filtered  noise  then  exponentially  modifies  the  nominal  target  strength 
which  is  our  output.  This  arrangement  allows  for  simulation  of  virtually  all 
types  of  scintillation  because  it  allows  control  of  both  the  statistics  of  the 
noise  and  the  characteristics  of  the  filter. 
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Pseudo -random  noise  is  generated  by  the  algorithm: 


Rl(n|-  65539  R1(n-»)<modulo23‘> 


(1) 


Rl(0»  * 12345 


Rn  -;{Ri(n)-<0-5lt231») 


(2) 


This  noise,  magnified  by  a scale  factor  K,  is  filtered  by  a digital  one -pole 
filter: 


X = C-X  . + C.KR 
n 2 n-  x 1 n 


where: 


At 


Cf  = 1 - e 


At 


Thus: 


At  A 

C,  = e * f~  t . 

C T 


x = ( 1 - ~)  ♦ X , + ~ • K • R 
n ' t ' n-1  t n 


The  values  selected  for  implementation  were 

t = 1.024  second 
At  =8  milliseconds 
At/T  = 2/256 
K = 128 


(3) 


(4) 
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Therefore: 


(5) 


x = HrX  4 + R 

n 256  n-1  n 

The  value  of  both  XR  and  the  eclipsing  constant  are  added  to  the  exponent  of 
the  normal  target  strength.  The  range  for  this  exponent  is  limited  to  0 through 
+63. 

3.  STATISTICS  OF  THE  ALGORITHM 
The  random  numbers  generated  in  equations  (1)  and  (2)  have  the  following 
uniform  distribution  with  mean  of  0 and  variance  equal  to  1/3. 


p(R  ) = 1/2  for  - 1 < R < + 1 
n n 


where  p(* ) is  the  probability  of  {• ). 


mean 


/•+* 

= PR  =J  R • P(R)  dR  = 0 


variance 


= aR 2 =/  (R  - PR)2  p(R)  dR  = 1/3 

CO 


After  passing  through  the  filter,  the  mean,  variance,  and  autocorrelation 
functions  of  any  one  sequence  [Xq.X^,  . . . XJ  as  derived  in  Appendix  B 
are  as  follows: 


t*X  = = 0 


(6) 


2 “V  2 k2ci  2 k2ci  < 

°X  ~ . ^2aR  ~ 2 + C . aR  “ 2 + C ‘ 3 


1 -C. 


1 


(7) 


1 
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E(X  , X . ) , 

_„  t n n-k  „ k 

P<k|  = a • a ~ * C2 

n n-k 


(8) 


Substituting  values: 


,x=° 


ax  = 21.4 


p(k)  = (0.  9921875) 


Time  = k • At  = 8 • k milliseconds 


(9) 


where 

k = number  of  samples 

The  amplitude  distribution  of  X is  Gaussian  with  the  above  mean  and  variance. 
These  perturbations  are  applied  to  the  exponent  of  target  strength  so  that,  in 
effect,  the  logarithm  of  strength  is  varied.  Applying  an  inverse  log  transfor- 
mation, we  find: 


1 2CTX 

p(X)  = e 


v/irr  i 


(10) 


p(S)  = 


v/irr  ♦ • S 


(InS  - InS) 
, 2 

1 2aX 

e 


for  -«o  < S < + » 
where  S = target  strength 

S = mean  of  S or  nominal  target  strength 


(ID 
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This  is  a log-normal  distribution  whose 

2 

CTX 

- R « 2 ^ o' 


«V  °R2 


1 - c. 


mean  = S e 


-2 

and  variance  = S 


= S e 


(12) 


(13) 


Note  that  the  average  target  strength  changes  from  its  nominal  value  S. 
Qualitatively,  this  can  be  understood  by  examining  a case  where  the  exponent 
of  target  strength  is  alternately  changed  by  +1  and  -1.  The  average  value  for 
the  exponent  for  this  case  is  unchanged.  But  when  one  inspects  the  corre- 
sponding strengths,  they  become  1-1/4  rather  than  1. 

4.  THEORETICAL  MODEL 

Radar  scintillation  discussions  invariably  refer  back  to  Swerling's 
classic  paper  [Ref.  2].  For  his  cases  I or  II,  which  apply  to  aircraft,  the 
dist-ibution  of  reflected  energy  is: 

_ S 

p(S,  S)  = | e S for  x ? 0 (14) 

S 

= 0 elsewhere 


where  S = input  signal-to-noise  ratio  (target  cross-section) 
S = average  of  S over  all  fluctuations 


8 


iwwywiniwMww 


1 v*  c-  *>'"rv'r  ■- y » ' v;  - - 

^-«^w’Wj.^“gagTygCTgSSig^ 


We  must  therefore  generate  a distribution  which,  when  applied  to  the  logarithm 
of  target  strength,  will  result  in  Eq.  (14).  If  p(x)  is  this  distribution,  then 


x = In  S or  S = eX 


dx  = 


dS 
S ~ 


e'X  dS 


P\x) 


E[S) 

dx 

dS 


S=e 


x 


S 

S 


x 


e 


p(x)  = 


x 

e for  - 00  < x < +oo 


(15) 


The  mean  and  variance  of  this  distribution  p(x)  as  derived  in  Appendix  C are: 


mean  = u.  =x-C  = lnS-C 
x 


where  C = Euler's  constant 

2 ir 

variance  = a = y 
x 6 


5.  COMPARISON  OF  MODELS 


A plot  of  the  theoretical  versus  the  implemented  distribution  is  included 
as  Figure  2.  Figure  3 illustrates  the  input  distribution  used  to  generate  the 
distribution  in  Figure  2,  Figure  4 illustrates  the  transform  used  to  convert 
theoretically  correct  noise  from  the  available  Gaussian  noise.  The  algorithm 
closely  realizes  theoretical  scintillation  for  valuer  greater  than  0.6. 


Figure  3.  Distributions  at  Input  to  Produce  Actual  or  Theoretically  Correct  Scintillation 


Figure  4.  Transformation  for  Samples  from  p(x)  to  Samples  from  p(y) 
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APPENDIX  A 

STATISTICS  OF  NOISE  GENERATOR 


The  basic  random  number  generator  is  the  sa:--e  one  used  in  the 
XDS  standard  FORTRAN  Numerical  Subroutine  Package.  It  generates  pseudo- 
random numbers  uniformly  distributed  between  1 and  positive  full  scale 

(X'7FFFFFFF').  Only  odd  numbers  are  generated.  Treating  the  generated 

31 

sequence  R^  as  numbers  scaled  over  2 


R 


1 


= 1 
= 0 


for  0 « Rj  < + 1 
elsewhere 


(Al-1) 


By  equation  (2),  these  are  converted  to  a uniform  distribution  in  the 
region  -1  to  +1: 


R 


for  -1  s R < + 1 
elsewhere 


(Al-2) 


Hence  the  mean  is: 


Hr 


J Rp(R) 

•CO 


. + 00 


dR  = 


/ i • R dR 


[< 


= 0 


(A  1-3) 


and  the  variance: 


r+® 

= 1 


ffR  = j-=o  P(RKR-HR) 


f r ^ i 

2,d  / JL  . rj 2,d  1 . R 1 

dR  “ J:  7 R dR  = I Lj"  J-l  = T 
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APPENDIX  B 


STATISTICS  OF  SCINTILLATION  SEQUENCES 


Each  sequence  of  scintillation  factors  is  generated  by  the  recursion 
relation  of  Eq.(3).  Expanding  this  equation: 

Xn  . = C . KRn  - 
0,  J 1 0,  J 

X,  . = C,KR , . + C.C .KR„ 

J 1 1,J  1 2 0,j 

X?  = C.KR,  + C.C,KR.  . + C.C,2KRn  . 

2,  j 1 2,  j 12  l,j  12  0,  j 

X,  . = C , KR x . +C,C,KR,  .+C,C,2KR.  . + C,C -?KRn  . 

3, j  1 3,  j 12  2,  j 12  i,j  12  0,  j 


X . = C.KR  . + C . C..KR  , . + 
n,j  1 n,j  12  n-l,j 


c.c"kr_ 

12  0,  j 


or 


X . = C.K  Y C"’1R.  , 
n,  J 1 L,  2 i,  j 


i=0 


.th 


The  mean  of  the  j sequence  is  obtained  by  summing  all  terms  in  the  sequence 


‘j'  sciK[(c 

n— oo  u 


|ij  - — i.,n  lie,0  + C2‘  + C2  + •••  CjjR^j 


(C“ 


+ c2*  + c\  + 


z)Ri  ,J 


Hi.  = i C.k[c?0  + c}  + •••  c"l  rRft  . + R.  . + R_  . + . . . R .1 
J n 1 L 2 2 2 J L 0,  J 1,J  2,  J n,  jJ 


n-*oo 


Preceding  page  blank 
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In  an  analogous  manner  we  can  find  the  expected  value  of  displaced 
samples : 

2 2 m 
C1KC2  2 

E<Xn  * Xn-m>  = 

1 ‘ ^2 


where  m is  the  displacement. 

3y  definition,  the  correlation  coefficient,  p,  is: 


p(m)  = 


E(X  • X ) E(X  • X ) 
1 n n-m  _ * n n-m 

(Tv  • (Tv 

X A OV 

n n-m  X 


which,  in  this  case  is; 

P (m)  = C™ 

The  mean  and  variance  of  he  transformed  sequences  may  be 

computed  fromEq.  (11),  winch  is  derived  in  the  following  manner.  Gaussian 

2 

pertubations  characterized  by  and  <r^ : 


p(X)  = 


(x-pxr 


are  transformed  by  an  inverse  log  transformation: 

S -•  eX  or  X = In  S 


where 


dS  X 1 
dX  ~ e ' S 


S = transformed  target  strength 
X - log  strengtn  = contents  of  target  files 


<x-'x> 


p(X) 

i 

i 

; dS  l 

,dx| 

X-lnS 

\Zv  o- 

X 

1ft 

Substituting  into  (A3-2)  and  (A3-3) 


_ y 

i r00  * 

•s.  = —J.  * e 

e 0 

, /■«  ~X  y 

(i2  = — y j e 6 (In  x)  dy 
e 0 

The  standard  forms  of  these  integrals  are  [from  Ref.  4]: 

I e”Px  In  x dx  = - — (C  + lnp) 

0 F 

r*-p*(  lnx)2dx  = CZ+4 
C 6 

where  C = Euler's  constant  = 0.577  

Using  these  and: 

it  can  be  shown  that: 

>*x  -hy**  "c+,n4> 

\e  / e 

- -(C  - In  e5^)  - x - C 

M ^ « MM  M ^ IT 

and  * ^ ’ x -**C  + C + ■£ 

8 (s-c)'  + r 


(A3-6) 

(A3-7) 

(3.  711,  R ;f.  3) 
(3.  714,  Ref.  3) 


(A3-8) 
(A  3-9) 

(A3- 10) 
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